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ABSTRACT 

In  this  paper  three  distinct  formulations  of  the  Pearson  chi-square  test 
statistic,  each  of  which  employs  random  interval  boundaries,  are  considered  in 
presenting  an  answer  to  the  following  problem:  Given  a  finite  sample  size  and 
specified  fixed  composite  null  and  alternative  hypotheses,  which  method  of 
estimation  of  the  unknown  parameters  leads  to  a  test  which  is  most  sensitive 
to  departures  from  the  null  hypothesis  in  the  direction  of  the  fixed 


alternative?  The  proposed  answer  is  based  on  a  finite  sample  size  approxima¬ 
tion  to  the  distribution  of  the  appropriate  quadratic  form  under  the  alterna¬ 
tive  hypothesis. N  In  this  investigation  the  estimators  considered  include  a 


simplified  least-squares  estimator  previously  developed  (Muhly  and  Gurland 
(1984),  MRC  Technical  Report  #2792)  and  others  satisfying  some  general 


regularity  conditions. 
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SIGNIFICANCE  AND  EXPLANATION 


There  are  many  tests  in  statistics  for  assessing  the  goodness  of  fit  of  a 
specified  family  of  distributions.  One  of  these,  the  Pearson  chi-square  test 
is  widely  used.  Different  formulations  of  this  test  are  possible,  with 
different  estimators  of  unknown  parameters  incorporated  in  it.  A  statistical 
question  of  primary  interest  is  the  following:  Suppose  a  finite  sample  is 
observed,  the  null  hypothesis  is  false,  and  some  fixed  alternative  hypothesis 
is  true.  Which  method  of  estimation  of  unknown  parameters  leads  to  a  test 
that  is  most  sensitive  in  detecting  the  falsity  of  the  null  hypothesis? 

Answers  to  this  question  are  particularly  important  when  testing  the 
reliability  of  a  product  where  some  underlying  family  of  distributions  is 
assumed  for  its  performance.  Another  important  situation  is  in  the  evaluation 
of  the  performance  of  some  treatment  when  an  alternative  treatment  is 
possible,  and  it  is  desired  to  detect  poor  or  superior  performance  of  the 
treatment,  as  the  case  may  be. 

In  the  present  paper  both  the  null  and  alternative  hypotheses  are 
composite,  and  the  investigation  of  the  above  question  is  based  on  the 
asymptotic  non-null  distribution  of  the  test  statistic  which  is  developed  here 
under  rather  general  conditions.  Estimators  considered  include  simplified 
least-squares  estimators  previously  developed  (Muhly  and  Gurland  (1984),  MRC 
Technical  Report  #2792),  and  others  satisfying  some  general  conditions.  The 
usefulness  of  these  results  is  in  assisting  the  statistician  to  make  decisions 
in  the  choice  of  the  formulation  of  the  Pearson  chi-square  statistic  and  of 
the  estimators  of  the  unknown  parameters  of  the  family  of  distributions  being 
tested . 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


EFFECT  OF  SOME  ESTIMATORS  ON  A  LARGE  SAMPLE  APPROXIMATION 
TO  THE  NON-NULL  DISTRIBUTION  OF  THE  PEARSON  CHI-SQUARE  STATISTIC 

Alan  E.  Muhly  and  John  Gurland 

1.  Introduction 

In  Chernoff  and  Lehmann  (1954)  the  limit  distribution  of  the  Pearson  chi-square  test 
statistic  is  obtained  when  unknown  parameters  are  estimated  by  the  method  of  maximum  like¬ 
lihood  based  on  the  full,  i.e.  ungrouped,  sample.  Since  the  full  sample  maximum  likelihood 
estimate  (m.l.e.)  is  generally  more  efficient  than  the  minimum  chi-square  estimate  (based 
on  the  grouped  sample),  the  use  of  the  full  sample  m.l.e.  should  lead  to  a  more  powerful 
test. 

Attempts  to  demonstrate  the  validity  of  this  conjecture  have  relied  on  some  notion  of 
the  asymptotic  relative  efficiency  of  the  Pearbon  chi-square  test  statistic.  In  Chibisov 
(1971)  and  again  in  Moore  and  Spruill  (1975)  it  is  shown  that  in  terms  of  Pitman  efficiency 
it  is  not  possible  to  assert  that  the  use  of  the  full  sample  m.l.e.  or  that  the  use  of  the 
minimum  chi-square  estimate  is  to  be  preferred.  Spruill  (1976),  using  the  notion  of 
approximate  Bahadur  slope,  also  showed  that  neither  method  of  estimation  is  superior  in  all 
situations. 

Although  these  studies  provide  many  useful  insights  into  the  behavior  of  the  Pearson 
chi-square  test  statistic,  they  leave  unanswered  the  question  of  primary  interest  to  the 
statistician.  In  practice  a  finite  sample  is  observed  and,  if  the  null  hypothesis  is 
false,  some  fixed  alternative  hypothesis  is  true.  Thus,  the  central  question  to  the 
statistician  is  this.  Given  the  finite  sample  sire,  if  the  null  hypothesis  is  false  and 
some  fixed  alternative  hypothesis  is  true,  which  method  of  estimation  leads  to  a  test  which 
is  more  sensitive  to  departures  from  the  null  hypothesis  in  the  direction  of  this  fixed 
alternative? 
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In  general,  if  the  null  hypothesis  is  false,  the  fixed  alternative  hypothesis  which  is 
actually  true  and  the  finite  sample  size  distribution  of  the  Pearson  chi-square  test 
statistic  under  this  alternative  hypothesis  are  both  unknown.  Thus,  it  is  not  possible  to 
give  a  precise  answer  to  the  question  raised  in  the  preceding  paragraph.  However,  the 
following  method  for  assessing  the  power  of  the  Pearson  chi-square  test  can  provide  an 
approximate  answer  to  this  question.  First,  suppose  that  a  random  sample  of  size  N  from 
a  population  with  unknown  distribution  function  F { • )  is  qiven  and  that  it  is  desired  to 
test  the  composite  null  hypothesis 

(1.1)  HQ  :  F<  •  )  =  Fq<*  ;0)  . 

In  (1.1)  it  is  assumed  that  the  functional  form  of  Fq  is  completely  specified  while  the 
parameter  vector  9,  which  is  an  element  of  a  set  0  contained  in  s-dimensional  Euclidean 
space,  is  unknown  and  unspecified.  Next,  if  Hq  is  not  true,  assume  that  the  fixed 
alternative  hypothesis 

(1.2)  H1  :  F(  •  )  =  F  (  •  ;  C  ) 

is  true  where  the  functional  form  of  F^  ia  completely  specified  and  the  parameter  vector 
5,  which  is  an  element  of  a  set  H  also  contained  in  s-dimensional  Euclidean  Bpace,  is 
unknown  and  unspecified.  Then,  since  the  Pearson  chi-square  test  statistic  may  be 
expressed  as  a  quadratic' form,  based  on  the  non-null  limit  distribution  of  the  random 
vector  appearing  in  this  quadratic  form  obtain  a  finite  sample  size  approximation  to  the 
distribution  of  the  test  statistic  under  Hj.  Finally,  from  this  approximation,  compute 
the  approximate  power  of  the  Pearson  chi-square  test  of  Hq  against  H^.  By  considering 
various  methods  for  estimating  0,  and  by  considering  a  number  of  different  specifications 
of  F,  in  (1.2),  it  is  possible  to  obtain  considerable  information  concerning  the  relative 
merits  of  competing  methods  of  estimation  in  the  test  of  Hq. 

In  this  paper  three  distinct  formulations  of  the  Pearson  chi-square  test  statistic, 
each  of  which  employs  random  interval  boundaries,  are  considered.  Based  on  the  procedure 
outlined  above,  finite  sample  size  approximations  to  the  non-null  distribution  of  each  of 
these  statistics  will  be  obtained  for  various  estimates  of  the  parameter  vector  9.  While 


the  test  statistics  considered  are  asymptot ical ly  equivalent  under  HQ/ 


it  will  be  seen 


that  this  is  not  necessarily  the  case  under  H-j  . 

Throughout  most  of  this  paper  it  will  be  assumed  that  both  Fq  and  .epend  only 

on  unknown  location  and  scale  parameters.  In  this  case  it  is  possible  to  define  a  least 
squares  estimate,  3^,  of  9  which  is  easy  to  compute  and  is  asymptotically  equivalent  to 
the  minimum  chi-square  estimate  under  Hq  (see  Muhly  and  Gurland  (1984a)).  In  view  of  the 

a, 

Chernoff  and  Lehmann  conjecture,  the  estimate  will  receive  specific  attention.  In 

addition,  other  estimates  of  9  are  considered  which  satisfy  fairly  general  regularity 
conditions.  In  particular,  Theorem  1  of  Dahiya  and  Gurland  (1973),  which  gives  a  finite 
sample  size  approximation  to  the  non-null  distribution  of  one  of  the  statistics  considered 
here  when  the  sample  mean  and  variance  are  employed  to  estimate  9,  is  obtained  as  a 
special  case  of  the  general  theory  developed  below  under  less  stringent  regularity 
condi t ions . 

The  next  section  summarizes  the  conventions  which  will  be  used  and  presents  the 
assumptions  that  will  be  required  in  this  paper.  In  section  3  the  three  formulations  of 
the  Pearson  chi-square  test  statistic  and  the  estimate  are  defined,  and,  in  section  4, 

some  preliminary  results  are  obtained.  Then,  in  section  5,  the  non-null  limit 
distributions  of  the  random  vectors  appearing  in  these  statistics  are  derived  and,  in 
section  6,  the  finite  sample  size  approximat ions  to  the  distributions  of  these  statistics 
under  H1  are  presented.  Finally,  in  section  7,  an  extension  is  discussed  which  permits 
consideration  of  some  parameters  not  of  the  location  and  scale  variety. 


2,  Conventions  and  Assumptions. 

The  following  conventions  will  he  adopted  in  this  paper.  First,  all  vectors  appearing 
here  are  column  vectors  and,  if  A  is  a  vector  (or  matrix).  A'  denotes  the  transpose 

V? 

of  A.  In  addition,  if  A  is  a  vector,  then  lAl  -  (A'A)  .  Second,  if  Hq  is  true. 


9^*  will  be  used  to  denote  the  true  value  of  9  and,  if  is  true  £  will  be  used  to 


denote  the  true  value  of  £.  Furthermore,  the  probability  measure  associated  with 
F  (  x  ;  9  ^  )  will  be  denoted  by  Pq  and  the  probability  measure  associated  with  F  (x;£  ) 


is  a  sequence  of  random  arrays  such  that  the 


will  he  denoted  by  Pj.  Third,  if  XN 

elements  of  f(N)XN  converge  to  zero  in  P^d-O  or  1)  probability  as  N  tends  to 

infinity,  the  notation  X  =*  o  ( 1/f  ( N )  )  will  be  employed.  Similarly,  if  XN  is  a 

t 

sequence  of  random  arrays  and  if  each  element  of  f(N)X^  remain  bounded  in  Pj 

probability  as  N  tends  to  infinity,  the  notation  X  =  0  (1/f(N>)  will  be  used. 

Jt 

Finally,  if  the  p  dimensional  random  vector  X  has  the  multivariate  normal  distribution 

with  mean  vector  u  and  covariance  matrix  t ,  this  will  be  denoted  by  X  ~  Np(u,^)>  In 

addition,  if  XN  is  a  sequence  of  random  vectors,  if  X  ~  Np(u,$),  and  if  XN  converges 

in  distribution  to  X  as  N  tends  to  infinity,  this  will  be  denoted  by  X  N  (u,J|). 

N  p 

In  order  to  obtain  the  results  given  in  sections  3  through  6,  it  will  be  necessary  to 
impose  regularity  conditions  both  on  Fq  and  on  F^.  The  basic  assumptions  which  must  be 
satisfied  by  Fg  are  the  following. 

Alg:  There  exists  a  parameterization  of  Fq  such  that 

1)  0  =  {0  =  {©  1 , 02  )  •  s  -*•  <  01  <  +“  and  02  >  0}  ,  and 

ii)  Fq  ( x ;  0 )  =  Fq(01  +  9  2  x )  where  FQ(x)  =  Fo(x>0)  for  0  =■  (0.1)*  . 

A2q :  Fq(x)  is  continuously  differentiable  in  x  and  fg(x)  "  ^  Fo'x)  I®  positive 

for 

all  finite  x. 

In  addition  to  these  assumptions  it  will,  at  times,  be  necessary  to  assume 

A3q :  J t4fQ(t)dt  <  “,  E0X  =  xdFo(x;0°)  -  _8°/02'  and  E0(X  +  91/92)2  = 

Ox  +  0°/0°)2dFo(x)00)  =  ( 1/9°  )2  . 

The  distribution  function  F ^  will  be  required  to  satisfy  similar  assumptions.  These  are 
given  by 

A 1 ^ :  There  exists  a  parameterization  of  F^(x;C)  such  that 

i)  =  »  {f,  =  ( ^  1  '^2  ) ' :  -"  <  ?1  <  ■*“”  an<1  5  2  >  '  and 

ii)  F  (x»£)  =  +?2X'  «bere  F^x)  «  F^xjS)  for  £  =  (0,1)’  , 

A2 y :  F^(x)  is  continuously  differentiable  in  x  and  f^(x)  =  Fj(x)  is  positive 

for  all  finite  x, 


and,  when  necessary, 


A3,:  J^,  t^ltjdt  <  »,  EtX  -  xdF^xiS  )  =  <]^2'  and  E1<X  +  “ 

Jlc  <*  +  tJ/c’^dP^xjC1)  =  ( 1/^2  >2  . 

In  addition  to  the  regularity  conditions  which  Fq  and  F1  must  satisfy,  it  will  be 

necessary  to  restrict  the  class  of  estimates  of  8®  as  well.  Basically,  in  order  to 

obtain  approximations  to  the  null  and  non-null  distributions  of  the  Pearson  chi-square  test 

statistic,  if  8  is  an  estimates  of  8^  it  is  only  necessary  to  require  that 
N 

y— -  *  0 

i)  /N  (9n"8  )  converges  to  a  proper  limit  distribution  under  HQ ,  and  that  ii)  there 
exists  a  5  in  H  such  that  /n  (8n-5)  converges  to  a  proper  limit  distribution  under 
H,.  The  approximations  given  below,  however,  also  require  that  these  limit  distributions 
be  normal.  Thus,  although  more  general  conditions  could  be  employed,  for  the  sake  of 
simplicity  the  conditions  which  estimates  of  9°  must  satisfy  are  given  in  the  following 
definitions . 

Definition  1 ;  An  estimate  8^  of  9°  is  said  to  satisfy  condition  Aq  if  there  exists  a 
function  h1^  :  R  *  such  that 

1)  h°  ( xldFg  (x;  9°  )  •  0  and  J  h° ( x  ) h° t x )  '  dFQ { x ; 8°  )  -  G°  where  the  elements 
of  G°  exist  and  are  finite  and  G°  Is  positive  semi-definite,  and 

«  0  1  N  0  - 

11)  9^-9  "  1  h  (*<,)  +  °p  'V^N)  where  X,,...,XN  denotes  the  original  sample. 

N  a-1  0 

—  *  0 
Definition  2i  If  5  i»  *n  element  of  z,  an  estimate  8n  of  8  is  said  to  satisfy 

—  *  2 
condition  A^C)  if  there  exists  a  function  h_  :  R  ♦  R  such  that 

i)  J ^  h_(x)dF  (x>5*)  ■  0  and  J  ^  h_(  x )h_(  x ) 'dF  ( x ;5 1 )  “  G_  where  the  elements  of 
5  5  5  5 

*  # 

G_  exist  and  are  finite  and  G_  is  positive  semi-definite,  and 

5  5 

li)  9*-5  -  i  l  h*(Xa)  ♦  Op  (1//N). 
a- 1  5  1 

_  * 

Definition  3:  If  £  and  £  are  two,  possibly  distinct,  elements  of  *•,  and  if  9^  and 

**  0  *  ** 

9^  are  two  estimates  of  9  ,  then  9^  and  9^  jointly  satisfy  condition 

-  *  2  **  2 

A„ (£ ,£ )  if  there  exist  functions  h  :  R  ♦  R  and  h  :  R  ♦  R 


- 


such  that 


li  >t  satisfies  condition  A  (F  )  with  ft  -  F  =  —  )  h(X)+o  (1//N)  and  G 

N  1  N  N  '•.-=-<*  P .  — 

0=1  F  1  F 


h_(  x  )h_(  x  )  '  il F 1  (  x  ,  F  ) , 


*  *  ~  ^  ^  *  — 
ii)  c'_,  satisfies  condition  Ai(F  )  with  0M  -  F  =  —  £  h  (X  )  +  o  (1//N)  and 

N  a=1  1  a  p- 


N 


1 


N 


1 


I  h  (x)h  ( x )  'dF  ( x ;  F, 1  ) ,  and 

J  — CO  -w  -v  1 


The  elements  of  the  matrix  G 


and  the  matrix 


G__  '  G 

l  l 


*  **  i 

h_(x)h^  (x)dF  <x;£  )  exist  and  are  finit 

f,  C 


is  positive  semi-definite. 


The  following  comments  are  in  order.  i)  The  parameterization  given  in  assumption  AIq 
and  A 1 ^  is  not  standard  for  distribution  functions  which  depend  on  location  and  scale 
parameters.  It  is  employed  here  for  the  reason  that  the  estimate  §N  given  the  next 

section  is  defined  in  terms  of  this  parameterization.  ii)  Since  the  problem  of  estimating 

0 

Q  is  secondary  to  the  problem  of  testing  HQ,  how  F0  is  parameterized  is  not  a 
significant  consideration.  Thus,  for  example,  any  distribution  function  of  the  form 

F(:x-u?  'd),  where  <  u  <  +«  and  o  >  0,  satisfies  A10  with  01  =  -p/a  and 

* 

a,  =  l/o.  Furthermore,  if  assumptions  A3q  and  A3<j  are  satisfied,  and  if  9^  is  given  by 
:  2 .  1  )  below,  thei.  it  is  shown  in  section  5  that  the  same  approximation  to  the  non-null 
distribution  of  the  Pearson  chi-sauare  test  statistic  is  obtained  as  when  Fq  is 
pirame* eri zed  in  terms  of  its  meand  and  variance  and  the  sample  mean  and  variance  are  used 
to  est  ijnat. e  these  parameters.  iii)  If  assumption  A3q  is  satisfied  then  a  natural  choice 
f f : r  an  estimate  of  is  aiven  by 


'2.1) 


\  -  1/^>’ 


■  1  N  2  1  N  -  2 

X  -  -  '  X  and  S  -7-  )  (X  -  X  )  .  Clearly,  9  satisfies  An  and,  if 

M  *1  *  a  N  N  '•  a  N  N  u 

■  t  "  1  HI  =.  1 


.1  n:mnt  t  i  on  Al,  ts  satisfied,  satisfies  condition  A  (f  ) 


However,  in  what  follows, 


with  the  exception  of  some  specific  examples#  it  is  not  required  that  A3q  and  A3 ^  be 

★ 

satisfied  nor  is  it  required  that  9^  be  given  by  (2.1). 


3.  Three  Forms  of  the  Pearson  Chi-Square  Test  Statistic  and  the  Estimate  3  • 

The  Pearson  chi-square  test  statistic  measures  deviations  from  Hp  by  considering  the 
difference  between  the  empirical  distribution  function,  FN(x),  and  the  (possibly 
estimated)  hypothesized  distribution  function,  F  (x;9°),  evaluated  at  K-1  points 
interior  to  the  support  of  Fq  .  Here  the  empirical  distribution  function  is  defined  as 


(3.1)  v*>  “if  *  W 

a=1 

where  ip  (y)  equals  one  if  y  is  less  than  or  equal  to  x  and  zero  otherwise. 

As  considered  here,  the  points  where  the  difference  between  f'N  and  Fq  is  observed 
are  obtained  as  follows.  First,  prior  to  observing  the  sample,  partition  the  interval 
[0,1]  into  K  (>  4)  intervals  and  let 

n.2,  o  -  6°0  < «“  <...< 1 

denote  the  endpoints  of  these  intervals.  Next,  for  i  =  1,2,..., K-1,  define  the  functions 


q1  (0  )  by 


g.<9>  =  F'hfiJjO)  , 


and  let  q  (0)  =  -«>  and  a  (0)  =  +°°  for  all  0  in  0.  Thus,  for  the  parameterization 
0  K 

-1  ,0.  „0 

specified  by  A1q,  g  ^  ( 6 )  =  for  i  =  1,2,...,K-1  where  =  FQ  (o^).  If  0 

were  known,  then  the  q.  (0°)  would  be  the  points  where  the  c  ifference  between  FN  and 

Fq  is  observed.  However,  since  0°  is  not  known,  obtain  a  preliminary  estimate  of  it, 

*  * 

9^,  and  compute  g^ ( 9^ )  for  i  =  1,2,..., K-1.  The  difference  between  FN  and  Fq  is 

* 

then  observed  at  the  K-1  random  points  g.(0  ). 

1  N 

For  0  and  0  in  0,  define  the  (K-1 ) -dimensional  random  vector  U(0,0)  by 

(3.4)  u.(0,0)  =  FN(gj(0>)  -  Fo(q.(0);0)  i  =  1,2 . K-1  , 

and  define  the  (K-1)  x  (K-1)  symmetric  matrix  r>  (0,0)  by 

(3.0)  D°('M)  =  ( F  ( g  .  .  <^T)  r9  )  [  1  -  F  (g  (^);9)!)K_1_  . 

'  0  ’mi n ( i , 7 ,  0  ’maxi i,))  1,3=' 


Note  that  if  0-0  then  FQ(q.(0);0)  =  6^,  which  is  known,  so  that 


-7- 


(3.6) 
and 

(3.7) 


U.  (9,9)  =  F  ( g .  (0)  )  -  6. 
1  N  l  1 


D°(0,0)  =  D°  »  (<5°  [1-6°  !  ) K” 1 

nun(i,j)  max(i,j)  a  r  j  =  1 


2  - 


Next,  for  9  and  6  in  0,  let  ^(9,9)  denote  the  quadratic  form 
(3.8) 


2  —  —  0  —  -1  - 
^(9,9)  =  NU  ( 9 , 9  ) '  D  (0,9)  U(9,9) 


If  denotes  the  preliminary  estimate  of  9  used  to  compute  the  points  g^(0N),  and 


*  *  0 

if  9  denotes  a  possibly  distinct  estimate  of  9  which  is  used  to  estimate 


F^(g.(°N^;^  then  it  is  shown  in  the  appendix  that  the  standard  definition  of  the 


2  *  *  * 

Pearson  chi-square  test  statistic  with  random  intervals  is  given  by  R  (9  ,9^  ). 


It  follows  from  (3.7)  that  the  true  value,  under  Hq,  of  the  matrix  appearing  in 

-1 


0  0  0-1  0 

(3.8)  is  given  by  D  (9  ,9  )  =  D  which  is  known  prior  to  observing  the  sample.  Thus, 


in  the  definition  of  the  Pearson  chi-square  test  statistic  it  is  not  necessary  to  estimate 


0 “ '  — 

D  by  D  (0^,0  )  .  Therefore,  for  9  and  9  in  0,  define  the  quadratic  form 


( 9  , 9  )  by 


-1 


(3.9) 


R^(0,9)  =  NU(0,9)'D0  U(0,9) 


*  **  o 

If  9  and  9  are  two  estimates  of  9  ,  an  alternative  definition  of  the  Pearson  chi- 
N  N 


•  2  *  ** 

square  test  statistic  is  given  by  RN^®N'®pj 


It  follows  from  (3.7)  that  when  =  ®N  '  =  Thus'  in  t*1*5  case 


the  finite  sample  size  approximations  to  the  distributions  of  these  statistics  are  the  same 


both  under  Hq  and  under  H-).  More  generally,  if  9  ^  -  9  “  o  (  1) ,  and 


if  9 


N 

*  •* 


satisfies  condition  Aq,  then  assumptions  Ain  and  A2n  are  sufficient  to  show  that 


IT(9  ,9  )  “  0  (  1//n)  (see,  e.g.,  Moore  and  Spruill  (op.cit.)).  Thus,  since  the  elements 

N  N  Pq 

of  D  (9,9)  1  arc  continuous  in  both  A  and  9,  RM(9  ,9  )  -  R  (9  ,9  )  =  o  (1). 

N  N  N  N  N  N  Pq 


However,  it  will  be  seen  in  section  6  that  the  finite  sample  size  approximations  to  the 


non-null  distributions  of  these  statistics  is  the  same  only  when  9..  *  9 

N  N 


Tn  order  to  define  the  third  form  of  the  Pearson  chi-square  test  statistic  considered 
here,  let  y  -  ( y , , . . . , yK_ , ) '  where 


(3.10) 


y.=  c.  +  I F  ( g . (9  ) )  -  S  ]/f„(c.  ) 
1  i  N  1  N  1  0  1 


and  0  satisfies  3  -  0  =  G_  (1//N).  Next,  for  0  in  0  define  the  ( K- 1 ) -d imens iona 1 

N  N  P 

0 


1 


'•„1 

■  .1 

1 


-H- 


random  vector  c(8)  by 


e ( 6 )  -  y  -  [1/9(6)19 
N 


where  1  is  the  ( K— 1 ) -dimensional  vector  of  ones  and  g(8)  =  ( g1 ( 0 gR_ 1 1 9  ))' .  In 
addition,  define  the  (K-1)  *  (K-1)  symmetric  matrices  n9  and  f  by 


n  =  diag[fQ (ct >, . . . ,fQ (cR_1 ) ] 

(3.12) 

r  «  ( n° )“ 1  d°  ( n°  > " 1  . 

*  2  * 

Finally,  for  9  in  0  define  the  quadratic  form  R  (8  ,8)  by 

N  N 

(3.13)  R^(8*,9)  =  NE(e)T_1e(8)  . 

**  0 

If  is  an  estimate  of  9  which  satisfies  condition  Aq ,  a  third  formulation  of  the 

*  2  *  *  * 

Pearson  chi-square  test  statistic  is  given  by  F.  (9  ,9  ). 

N  N  N 

It  is  shown  in  Muhly  and  Gurland  (op.cit.)  that  »6.,  )  =  0  (1) 

N  N  N  N  N  N  F 

*  Q  _  ** 

when  8  -  8  ■=  0  (1/»N)  and  8^  satisfies  Aq.  However,  it  will  be  seen  in  section  6 
0 

that  the  finite  sample  size  approximations  to  the  non-null  distributions  of  these  two 

*  ** 

statistics  are  the  same  only  when  8^  -  8^  =  o  (1). 

0  ~  1 

The  minimum  chi-square  estimate  of  8  ,  8^,  is  defined  to  be  any  estimate  for  which 
there  exists  a  sequence  of  positive  constants  p  ,  with  lim  =  0,  such  that 


(3.14) 


rn(vv  <  S2  +  p» 


where  8^  is  the  preliminary  estimate  of  6  used  to  compute  g^fS^)  for  i  =  1,2,..., K-1. 

•  n  —  ~  “  2  * 

If  0N  -  8  *  Op  (  1  // N )  ,  a  least  squares  approximation  to  0^,  which  minimizes  ^(9^,6) 

0 

with  respect  to  9,  is  given  by 

(3.15)  8  =  t  ( 1 ,  g<  e*)  )  •  r-1  ( 1  ,g(  8*  > )  r1  (  1  ,qr  ( 9* ) )  •  r“  V  . 

N  N  N  N 

If  3,  is  any  estimate  of  8°  satisfying  (3.14),  and  if  8,,  -  9^  =  On  (1//N),  then  (see 

N  -P0 

~  «  r-  2,*~  2*? 

Muhly  and  Gurland  (op.  cit.))  ”  %  =  °p  and  KN*9h,l3N*  ”  !<N^9N'^N*  ~  °p0  *  '  *  ' 


4.  Preliminary  Results* 

In  this  section  several  consequences  of  assumptions  A1^  and  A2^(Jl  -  0  r  1)  are 
obtained  which  will  be  used  in  the  remaining  sections.  First,  note  that  for  any  \  anc 


and  the  result  follows  from  corollary  5.1. 

The  non-null  limit  distribution  of  can  be  obtained  from  lemma  5.2  and 

N  N 

theorem  5.1.  However,  the  following  direct  argument  provides  a  more  convenient  expression 
for  this  result. 

*0  — 
Theorem  5,3:  Let  8n  denote  an  estimate  of  8  and  suppose  that  there  exists  a  £  in 

*  —  « 

5  such  that  0  satisfies  condition  A.(^).  Define  £  as  in  lemma  5.2  and 
N  1 

1  et 

m  (T,c>  =  i  -  n°(r,“)(i,g(£))E(?)(n0)'1  , 

MgK1,?,?)  -  -[M7(^,5)li1(C,£1)  +  n°(?J)]  (  1,g(£))  , 

and 

t(4)(f.1,T,!>  =  (m7(T,0  ;  MB(e1,?,5)]s\c1.r>  [M7(C,U  I  m8(C  1  ) ] •  • 

Tf  8^.  is  defined  by  (3.15),  and  if  Al^  and  A2^  are  satisfied  for  i  =  0 ,  1 
then ,  under  H ^  , 

^(U(VV  '  Bm(f1'T^)!  *  Vi(0't(4),?1^11 
1  ■—  50 

where  fl  ^  ,  h,h)  is  defined  in  theorem  5.1. 


-23- 


a  (  Z) 


2 

-  u)2  < c  )to  ^  (  Z  )  ]  /  [n>1ul  (c  )  -  (  c )  ]  , 

b(Z1  -  -  '02(c  )u>2  (  2 )  l/taiya  (c>  -  UJ  2  (  C  )  2  ] 

n i:lc  •  j  '  -ill  at  :  nn.i  show  that 

E(f.)Z  (  a  (  Z  )  +  b(Z)C1  b  (  Z )  £  ^  >  * 

*  0  - 
. :  Si-  *va»  9  is  an  estimate  of  9  for  which  there  exists  a  £  in  H 

—  N 

*  -w 

si:t  h  *■ : .  a  r  satisfies  condition  A  ^  ( £)  .  In  the  notation  of  corollary  5.2 

i  1  — 

i  n  f  , r  > 

•  2  > 

:  -  r  a  f  r  )  ♦  b  ( p  )  £  ^  ,  b  ( p  )  C  .  ) 1  • 

-  0-1 
f :  )  =  F<  f  )  i n  ) 

-  ((Ub(p))i2  -  E(r)(n°)'1n1(c,C1)(i#g(^’)  )3  , 

arid 

[m5(£)  •  M6tc 1  ,r>  i  •  . 

If  b(o)  >  0,  and  if  A1 1  and  A2 i  are  satisfied  then,  under  H^» 

oueN  -  a  *  n2 <o,t(3)<(\o? 

where  3  is  defined  in  (3.15). 

N 

It  follows  from  (3.10),  (3.15),  (5.31),  (5.32),  (5.35),  and  the  definition  of  p 
that 

%  -  E(%)y  -  e(9*hc  *  (n°)_1u(9*,0*)] 

-  7  +  F  (T)P  +  [E(0*)p  -  E(T)p)  +  (8-0 

N  N 

<■  E(0*)(n°)'1[U(0*,0*)  -  B(2)<*i\?>) 

=  ’  *■  (1+b(p))(0*-O  +  E(0* ')(H°  )-1  [U(8*,9*)  -  B.-.tfVo)  • 

N  N  N  N  \  Z  ) 


r  ur  the  r*mnre ,  (3.6),  (5.13),  and  (5.15)  show  that 


[i;  m4(51,T,5)  i  -M,g(5>]|[(FN<g<5>>  -  <5 1  (5  ,€ 1 ) )  *  •  (9*  -  5)'  ;  (6**  -  ?)■]•  +  op  (1//n) 

and  the  result  follows  from  lemma  5.1. 

1 

** 

It  is  interesting  to  note  that  (5.29)  shows  that  e(6N  )  can  be  expressed  as 

c(eN  )  -  y  -  F01(«°(9n,9n  )) 

0-1  *  0  *  **  0  *  **  0 

-  c  +  ( n  )  EFN(g(9N))  -  6  (9n,9n  )  +  6  (9N.9N  )  -  «  1 

‘  tF“1(«°(9*,e**))  ~  Fg\s°<5,5>)  +  F"1 (6° (T.e ) ) 

-  c  -  [  1 , g(T)  1 5  +  <n  )  U(9n,0n  )  +  (nU)  [6  (5,0  -  6U] 

*  [(n0)_1-n°(c,o_1]  [«°(6*,e**)  -60(T,oi  +  op  <i//n>  . 


-  -N.  *  ** 

Thus,  if,  in  theorem  5.2,  5*5,  i.e.  if  9-9  =  o_  ( 1 ) ,  then 


(5.30) 


e(9**>  -  (n0)‘1U(9*,e**)  +  op  (1//N)  . 


Therefore,  the  following  corollary  to  theorem  5.1  has  been  proved. 

*  **  o  *  *  * 

Corollary  5.4;  If  9^  and  9^  are  two  estimates  of  0  ,  if  eN  “  =  °p  ( 1) ,  if 

_  h  It  h  ^ 

there  exists  a  5  in  =  such  that  9^  and  9^  jointly  satisfy 
condition  and  if  Alj  and  A2 ^  are  satisfied  for  1=0,1 


then,  under  Hj, 

/n  (e  (0**  )  -  (n°)"1B<1)(51,5,5)]  $  ( 0,  (IT0)'1*  ( 1 )  (5  1 ,5,5  )  <n°  _1)  . 

It  remains  to  determine  the  non-null  limit  distributions  of  8„  and  U(9  ,8  ).  For 

N  N  N 

this  purpose,  define  the  2  x  ( K— 1 )  matrix  E(5)  by 


(5.31) 


E(5)  -  [d,g(5))  ,r"1(l,g(5))]'1d,g(5)),r'1 


where  5  is  an  element  of  H,  and  notice  that 


E(5)c  =  e(5)  [  1  ,g< 5)15  =  5 


Next,  for  Z  »  ( Z^ , . . . ,ZR_ 1 ) '  define 


=  iV’l  , 

u2(z)  =  1 1r-1z  , 


dlJ(Z)  =  cV’z  , 


P  (f. 1  ,r.  .1 )  =  c  -  [1,gKm  ♦  (H0)-1  ['  1  (£  ,£ 1  )  -  5°]  , 

=  [d2/T2)l  -  li1;?,?,1)!  M,g(S)]  , 

and  let 

t  {e\r,?>  »  [i  m4(c1,t>?)  j  -( i,q<r>  > i s 1  [i  |  M4(?1,f,i)  j  -  ( n  • 

If  A 1  ^  and  A2jj  are  satisfied  for  t  =  0,1  then,  under  H1# 

/n  ( c ( e** )  -  8(3)(f.1,r,U)d*  Vi(°,*(3)<e1,^))  • 

Proof ;  For  5  and  i  in  Z,  let 

f"1  <60<r.^)  )  *=  <f"1(S°U,S)) . Fo1(SK-1(^n' 

=  [l,g(T))C 

by  (5.3).  Then  it  follows  from  (3.10)  and  (3.11)  that 
e<9**>  «  y  -  (l,g(9*)ie** 

-  c  +  (n°)_1[FN(g(9*))  -  «°]  -  F~1(6°(9*,9**))  . 

Since  6°  ( 9*,  9** )  -  6°(?,|)  =  0D  (1//N)  (see  (5.14),  (5.16),  and  (5.17)),  and 
N  N  P  ^ 

since  —  F„ 1 C t )  l  „  =  1/*^(£.S),  lemma  4.2  shows  that 

dt  0  I  .Or'  l 

t«=o1(5,5) 

(5.20)  f”1<5°(0*,6**  )  )  =  F~1(5°(I,?))  +  n°(?,S)-1  [«°(9*,9**)  -  «°(C,I)]  +  op  (1//N)  . 
Furthermore,  it  follows  from  (5.13)  and  (5.15)  that 

FN(g(9*))  -  6°  =  [«  1  (TrC1 )  -  «°)  +  (FN(g<r>>  -  «°  (T,C  1 )  J 

-  nVT.c’Hl.gtTjnS*  -  ?)  +  op  (1//N)  , 

and  (5.14),  (5.16),  and  (5.17)  show  that 

«0«e:.O  -  «°(C,5)  =  n0(T,c)[i,g(?)]  (9**  -  l) 

N  N  N 

-  n°(T,5 )  [i ,g(T)]  (0*  -  T)  +  o  d//N)  . 

”  i 

Thus,  r(0  )  =  6  (C  ,  E,£)  + 

N  l  3  ) 


which  {5.20)  and  (5.21)  were  obtained.  In  general,  it  follows  from  Muhly  and  Gurland 

#  *  * 

(1984b)  that  if  and  0^  are  invariant  with  respect  to  the  group  of  linear  transfor¬ 
mations  with  positive  slope,  then  6  ^  ^  ,£  )  and  of  theory..  5.1  do  not 

depend  on  unknown  parameters. 


Two  further  applications  of  corollary  5.2  yield  the  non-null  limit  distributions  of 


e  (8*  )  and  e (3  ). 

N  N 


To  see  this  let  c  =  (c-j , . . .  ,cK-1  >  and  notice  that  for  all  5  in  = 


(5.25) 


=  c 


Then  it  follows  from  (3.10),  (3.11),  (5.8),  and  (5.25)  that 


(5.26) 


e(Sj  =  y-c  =  ©*)  . 

N  N  N 


Next,  let  P  denote  the  (K-1)  x  (K-1)  matrix 


(5.27) 


P  =  I  -  (1,g(0*))t(1,g(9*)),r‘1(1,g(0*)))“1(1,g(e*))T"1 
N  N  N  N 


It  is  shown  in  Muhly  and  Gurland  (1984a)  that  P  does  not  depend  on  9„  and,  since 

N 


P ( 1 » g ( ©  )  )  =  0,  Pc  =  0.  Furthermore,  it  follows  from  (3.11),  (3.15),  and  (5.26)  that 


(5.28) 


e(3N)  =  Py  -  P( n° ) _1O(9*,0* )  . 


Thus,  the  following  corollary  has  be-:n  proved. 

*  0  - 
Corollary  5.3:  Let  @N  denote  an  estimate  of  8  and  suppose  that  there  exists  a  5  in 

*  _ 

=  such  that  0N  satisfies  condition  A^S).  If  Al.,  and  A21  are  satisfied 
and  if  9^  is  defined  by  (3.15),  then  under  H^, 

/n  (e(0*)  -  <n° )_1e (2 } (C 1  ,T) )  -  NK_1(o,(n°)'1{(2)(51,T)(n0)'1) 


/n  ( e ( 9  )  -  potVb 


N>  -  p(Hu)"1b(2)(?,1,£))  *  NK_1(o,p(nu)'1i(2)(51,E)(n0)"1p- > 


For  general  0  the  followina  theorem  gives  the  non-null  limit  distribution  of 
N 

*  * 

e,9n  >• 

*  **  o 

Theorem  5.2:  Let  9_  and  0_  denote  two  estimates  of  0  and  suppose  that  there  exist 
'  N  N 

_  ^  *  *  * 

£  and  £  in  S  such  that  and  0N  jointly  satisfy  condition 


A.  <F.,£).  Let 


92(2  =  J-»<t  “I)  , 


then  ( '>.-»)  shows  that  the  ith  element  of  *s  9*ven  by 

6i  “  5i  1  “  , 

and  straightforward  calculations  show  that  the  (i,j)th  element  of  *s  9iven 


«\  ..1  +  [a’fc.)  +  a’(c  )]f  (c  ) 

max(i,j)  1  l  2  2  l  1  j 


(5.21)  +  [  a  (  c  .  )  +  —  a_(c. ))f ,(c. ) 

1  1  2  2  ]  1  l 


i.j  =  1,2,...,K-1 


+  :i  *  (Ci+Cj,q2,l/2  +  CiCjg22/4,f1<Ci)f1<Cj) 

In  theorem  1  of  Dahiya  and  Gurland  (op.  cit.)  the  following  variation  of  corollary  5.2 

i  55  obtained.  First,  it  is  assumed  that  Fq(x;9)  and  F^(x;£)  are  parametrized  so  that 
<5.2  2.)  Kr(x;fl)  =  Fq(  (x-0  1 )  //  6  2 )  ,  j ^  xdF0<x;9°)  -  0°,  and  j  ^  ( x-0® )  2dFQ  ( x  >9  °)  =  0°  , 


‘  ^  .  23  )  F,  (x;f)  =  F1  ( ( x-£  ,  )//iTJ) ,  jZ*.  xdF-|  (x)?1 >  -  e],  and  (x-5  ^ )  ^dF  1  ( x  ;£  ' )  =  £  j  • 

If.  addition,  it  is  assumed  that  j”!,  t^itldt  <  •  for  J.  =  0,1  and  that  0N  *  'xn'Sn''* 

.  * 

Then,  considering  the  K-dimensional  random  vector  defined  by 


VVV 


if  i  =■  1 


VCO  -  VVV  -  Vl'VO  if  i'2'3 . K-1 


-0  , ( 0* ,0* ) 

K-1  N'  N 


if  i  =  K  , 


and  under  more  stringent  regularity  conditions  than  required  here,  they  obtain  expressions 

*  * 

for  the  asymptotic  mean  vector  and  covariance  matrix  of  ^®n'9N^  “hich  are  easily  shown 
to  be  equivalent  to  (5.20)  and  (5.21).  Thus,  in  this  case,  how  F0  and  F1  are  para¬ 
metrize)  affect n  neither  the  null  nor  the  non-null  finite  sample  size  approximations  to  the 

list  r  >  hut  ion  of  '-'(9 *,"*)  . 

N  N 

Ar  other  fw.int  to  notice  about  (5.20)  and  (5.21)  is  this.  Neither  expression  depends 
on  ,wr  pa  r.imet  ,.r  n .  However,  this  result  is  not  restricted  to  the  conditions  under 


denote  the  distribution  function  of  N(0,Y'  »E,E)Y>  random  variable.  Then 

F^(-“)  ”  0,  F^(+*»)  •  1,  and,  by  theorem  5.1,  f^n(x)  converges  weakly  to  F^tx). 

-  *  *#  1  —  —  _ 

Furthermore,  sirce  /n  y' [U(9„,9  )  -  6,,.<E  ,E,E>1  <  2/NOC-1 )  y  <  »  for  all  N, 

n  n  U  ) 

F  „(-*•)  ”  0  and  F  (+»)  »  1.  Therefore,  by  theorem  8.6.2  (if  y’  t,,,(C,E\EJ T  >  0)  o 

T"  IN  (  1  ) 

by  theorem  8.6.1  (if  y'  Y  »  0)  in  Ash  (1972),  F^fx)  converges  to  F^fx) 

uniformly  in  x  as  N  tends  to  infinity.  Since  Y  was  arbitrary,  it  follows  that  for 

«  ** 

large  N  the  distribution  of  »N  U(8  ,9  )  can  be  approximated  by  the 

N  N 

Nr_i(/n  B(1,(51,I,C),  t(1)(?1,T.O)  distribution. 

*  ** 

The  following  corollary  to  theorem  5.1  treats  the  case  where  9^  =*  0^  . 

*  0  — 
Corollary  5.2:  Let  6  be  an  estimate  of  9  and  suppose  that  there  exists  a  £  in 
N 

*  _ 

such  that  9^  satisfies  A((E).  Let 

b(2)(e’,T)  «  s’d.s1)  -  6° 
m3(c1,e)  -  -n’d.e’Hi.gd)] 

and  let 

$(2)<c\T>  -  [i  ;  m3(e1,e)]s1(e1,ehi  :  m3(e1,e")]'  . 

If  A1 j  and  A21  are  satisfied  then,  under  H^, 

✓N<u<e*,e*)  -  s(2) (E1,?) )  ♦  4(2)(51,5))  • 

As  an  application  of  corollary  5.2  suppose  that  A3q  and  A31  are  satisfied  and  that 

*  *  i  « 

9^  is  defined  by  (2.1).  Then  9^  satisfies  condition  A^E  )  with  h  ^(x)  given  by 


(5.19) 


h  ,  ( x ) 

F  * 


Furthermore,  if 


(5.20) 


(l 

*  +  *y2 

(xH]/e’)2  -  cve’j2 

a J (x)  «  J  ^  tf f ( t  )dt  , 
a2(x)  -  Jl,(t2-1  )f 3 < t )dt 

^2,1  “  il.  and 


.»  2  2 
j_(t  -i  >  f )(t)dt 


and  let 


M2(5,5)  =  -n°(5.5)t1.g(5>)  , 

$(1)<c\r,?)  =  (i  i  M1<c1,r,l)  •  M2<r,e)]s1(c1,r,c)(i  •  j  m2 tr,c > j •  . 

If  A 1 ^  and  A2^  are  satisfied  for  1=0,1  then,  under  H^> 

/N(u(9*,e**)  -  6(1)(e1,r,l))  ♦  NK_1(o,t(1)(e1,|‘,C))  • 

Proof:  U*9N'9N  *  =  F  (g(0N> )  "  <$°(0N,9N  '*  Appealing  to  corollary  4.1, 

(5.13)  fN(<?(6N)>  =  FN(g(^n  "  +  51(®n-51)  +  Op^(1//N)  , 

and  corollary  4.3  shows  that 

(5.14)  6°(9N'9N*)  =  6°<9n'^)  +  ^  °  (^"'9  )  "  «°(5,1)  +  op  (1//N)  . 

Thus,  u<e;,9;*)  -  b(1)(£1,t,c)  -  (FN<g<?))  - « 1  (T,^1  > j  +  -s’tf.e1)] 

-  («°(0*,5)  -  «°(T,I>]  -  [5°  (T,0** )  -  «°(T,C)]  +  Op  (1/h). 

Next,  lemma  4.1  shows  that 

(5.15)  61  <e* .e1 )  -  « 1  (T. c 1 )  =  -n1(c',C1)ti.g(C)](0*-C)  +  op  (1/iT)  , 

and  that 

(5.16)  «°(9*,1)  -  5°(5.b  *  -n°<r.|)  [i.g(O)0*-O  +  op  ( 1//N)  . 

Furthermore,  corollary  4.2  yields 

(5.17)  6°(?,0**)  "  6°<^>  ”  n°(5,?)[1,g(?)]  (0^*-l>  +  op  (1//N)  . 

Thus,  U'9n,9n  *  ”  9(ij  ',£,?)  “ 

(I  :  M.(E\?,I)  •  H,(E,()][(FI<I(I)]  -  «  1  (F,£ 1 ) )  ’  - 

(5.18) 

<e*-T)  *  ;  (9**-C  )  *  1  *  +  op  (1//N) 

and  the  result  follows  from  lemma  5.1.  ■ 

In  theorem  5.1  it  should  be  noted  that  the  convergence  of  the  distribution  function  of 

/n(U(0„,0„  )  -  B,  , .  (51  ,T,X  )]  to  the  distribution  function  of  a  N  .  (0,$  . .  (5  1 ,£  ,%  ) ) 

N  N  (  1  )  1  V  •  ) 

random  vector  is  uniform  as  N  tends  to  infinity.  To  see  this  let  Y  be  an  element  of 

R*'1  subject  to  the  restriction  that  I Y *  <  ”  but  otherwise  arbitrary,  let  F^tx) 

*  •«  1  —  - 

denote  the  distribution  function  of  iN  l,'[U(®jI,®|)  )  “  8^j(t  ,£,£)],  and  let  F^(x) 
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and  let 


a  ’  1  v ,  for  r  and  f  in  -f  define  tne  ( K-1 ) -dimensional  vectors  6'(f,  ,f  )  (fc*  0,1)  ar 
fK-1)  *  f  K—  1 )  symmetric  matrix  P^(f,,£)  by 

.5'.  =  <**<r,(:> . *  -  o,i  , 


'5.6)  n  V,6>  -  (■*.  V  ,(6,0!  1  -  61  .<6.6)1 )?~1  .  . 

min(i,j)  max(ifj)  i ,  j  *=  1 

*  **  Q 

Lemma  5 . 1 :  Let  6n  and  denote  two  estimates  of  0  and  suppose  that  there  exist 

-  *  **  _  -N. 

£  and  f,  in  •:  such  that  0^  and  0N  jointly  satisfy  condition  A  (£  ,£  )  • 

*  ** 

Define  the  (K-1)  *  2  matrices  W_  and  W  by 

l  l 

W_  -  ,  _  (x)h_(x)’dF  ( x  j  c  1 ) 

6  g<6  )  5 

W  =  J°°  T  (x)h  (x)’dF  (x>6  ) 

6  q(K)  l 


1—1  * 

d  (6,6  )  w  w 


1  1  —  ~  • 

s  (6  ,6,6)  «  W 


G_  *  G 

6,6  6 


Then,  under  H1 , 


/Nt(F(g<6))  -  -51(6,51))*  :  ( e*  -  6  )  •  <9**  -  6)’  1  ’  $  N_,(0,s’(6  ’,6,6  ))  . 

N  *  N  N  K+j 

Proof:  Let  Z(x)  =  ! ( f  lx)  -  6 1 (^,6 1 ) ) ’  :  h  (x)’  h**(x)’J'.  Then 

g(6)  66 

mN(q(6))  -  F’fOF1  >)'  <9*-6)'  (9**-6)')'  *  i  l  ZOO  +  0p  (1//N).  Further 

a*  1  1 

more,  (i)  r  Z(x  )dF  ( x;6  1 )  =  0,  (ii)  Z(  x  )Z  ( x )  1  dF ,  (  x ;  6  1 )  =  S1(61,6,6),  and 

(iii)  if  a  /  P ,  Z(X  )  and  Z(Xa)  are  independent.  Thus,  by  the  multivariate 

01  p 

,  N  1  1  _  ~ 

version  of  the  Li ndeberq-Levy  theorem,  —  )  Z(X  )  ♦  N  (0,S  ) )  . 

/N  a=  1  a  K  3  . 


irol  1  a  ry  *> .  1  :  Let.  A  denote  an  estimate  of  A  and  suppose  that  there  exists  a  f  in 


su'-h  that  A.,  satisfies  condition  A*(£).  Define  W 


as  in  lemma  5.  1 


sup  lv(u)  -  <  Ej/2 C.  Let  denote  the  event 

|u-I.|<a 

*  *,  . _ 

If  Ajj  occurs,  then  g*  N  is  in  S(c^)  and 

'VVViV'^y  -  V*i*2*i«”"'*i«”'  ' 

*  _ 

lv(FQ(gi  N;C))  -  vtS^JI  <  e2/2C.  Thus,  the  occurrence  of  A^  implies  that 
Ci,N(5;T,|)  -  /N|f0(n,g;<N)?)(i,g;  N)(C-?)  +  r^.T.%) 

-  f0((i.gi(T))5)(i,gi(T))(C-C)  -  r(gi(I),I,?)| 

<  /NI5-CI  •  ‘f0(<1'9i(N)5)(1,g*>N)  -  f0((l,gi(?))i)d,gi(?))l  +  2  ✓n’  y  <15-1*  » . 

Since  y { x )  »  o(T)  as  T  ♦  0,  /N  y(C//N)/C  +0  as  N  ♦  «■  and  thus 

2/n  y(C//n)  <  e2/2  for  N  sufficiently  large.  Therefore,  for  N  sufficiently 

large,  the  occurrence  of  A^  implies  that  sup  N<C;€,?)  < 

K-£I<C//n 

/N(C/rN)(e2/2C)  +  2/N  y(C//N)  <  Ej.  Finally,  since  F0(g1Nie)  -  S'.  -  op  ( 1 ) , 
P^Ajj)  ♦  1  as  N  ♦  «. 

Corollary  4.3:  If  £  and  £  are  any  two  elements  of  if  9^  -  5  =  op  ( 1 )  and 

8**  -  l  -  Op  (1//n),  and  if  Al0  and  A20  are  satisfied,  then 

vve>6r>  ■  +  F0(9i(5)»5)  -  o  . 

5.  The  Non-Null  Limit  Distributions  of  U(0__,0._  ),  e(0._  ),  and  S._. 

"  N  N  N  "  ■  ■“  N 

All  limit  distributions  obtained  in  this,  and  the  next,  section  are  based  on  lemma  5.1 
and  corollary  5.1.  Before  stating  this  lemma  and  itB  corollary,  however,  it  is  convenient 
to  introduce  the  following  notation.  For  any  vector  a  »  (a1 , . . . ,aK-1 ) ' ,  let 

(5.1)  Y  (x)  -  (*  (x), ...,*  (x))'  , 

a  a1  aK-1 

and  for  any  £  in  3,  let 

1  « 

(5.2)  F„(g(5))  -  (FN(g1(5)),...,FN(gK_1(5)))'  =  -  \  Y  _  (X^)  . 

a=1  g(S) 

Furthermore,  if  ?  and  £  are  any  two  elements  of  3,  define 

(5.3)  «J(T,S>  -  Ft(gi(5)»5)  i  -  1,...,K-1,  i  -  0,1  , 

and  notice  that  when  £  ■  £,  (5.3)  becomes 

«J(C,i)  -  =  Fe<Cl > 
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(5.4) 


bounded,  on  B(E),  v(u)  is  continuous  and  bounded  for  u  in  B(E). 


Choose  £1  and  r 1  so  that  0  <  e1  <  8  and  0  <  T 1  <  $2<  By  the  Cauchy- 
Schwa  rz  i nequa 1  it  y 

|r(x,8,9) 1/10-01  <  (1+x  )  2  |f0((1,x)9)  -  f Q( ( 1 ,x>6 ) | 
for  all  (x,8)  in  x  b(3,Xj).  Moreover,  since  A ( E  ^ ,x  ^  )  is  compact  and 

fg  is  continuous  on  Afe^x^,  fg  is  uniformly  continuous  on  A(e^,T^).  Thus, 
given  >  0  (with  <  cl'  there  exists  a  t2  >  0  (with  x2  c  x^)  such  that 

|(1,x)9  -  (1,x)S|  <  x  2  implies  |fQ((1,x)0)  -  f  Q(  (  1  ,x  )0  )  |  < 

Furthermore,  since  |<1,x)9  -  ( 1 ,  x  )3  |  -  a|  ( 1  ,x)  (0-3)  j  <  (  1+x2)1^2  10-01  <  M(e  )I0-§I 

for  all  x  in  D ( e  ^ ) ,  10-91  <  implies  that  |f  ((1,x)0)  - 

fQ(  (1,x)3)|  <  E2/M(e.()  for  all  x  in  D(e1).  Therefore,  if  19-91  <  t  2/m  (c^), 

sup  |r(x,0,0)  |/I0-9I  <  since  )  c  DIe^.  Thus,  sup  |r(x,9,9)|  » 

xesfe^  xesfe.) 

o(  1 9-3 1  )  as  1 9-3 1  +0.  „ 


Corollary  4.2:  Suppose  that  5  and  5  are  any  two  elements  of  S  and  let 


Fq (gi (£);£) .  If  A 1 g  and  A2g  are  satisfied,  and  if 


8  -5=0  (1//N),  then 

N  ♦*  1  •* 

Fo(gi(T);0N  )  =  J .  +  f0(€1+52gi(5))(1,gi(r))(8**-!)  +  op  (1//N). 


Lemma  4.5:  Suppose  that  5  and  5  are  any  two  elements  of  S  and  let 

<5 ^  “  F()(gi(5))5).  Let  N  be  a  sequence  of  random  variables  such  that 


F_ ( g ,  ;5)  -  i,  =  o  (1)  and,  for  5  in  E,  define 

u  i,n  1  f j 


?1/N(5;5,C)  -  ^|F0(gi<Nl5)  -  F0(gi>N;5)  - 

If  A10  and  A2q  are  satisfied,  then  for  any  C  >  0, 


sup 


5i/N(5,C,C)  =oPi(l) 


I5-5I<c//n 
ists  an  £1 

B(£^).  Furthermore,  lemma  4.4  shows  that  y(x)  =  sup 


Proof :  By  lemma  4.4  there  exists  an  £^  >  0  such  that  v(u)  is  bounded  and  continuous  on 


sup 


I  r  (x,5  ,5  )  |  “  o(x  ) 


I5-£|<x  xesi(e) 


as  x  ♦  0.  Choose  E„  >  0  and  find  an  a  <  e.  such  that 


Corol  lary  4.1:  Suppose  that  F^(x;r,  )  is  continuous  in  x  and  that  £  is  an  element  of 

Set  0  .  =  F,  la.  (f  )  '  )  and  let  g-  M  be  a  sequence  of  random 
111 

*  1  — 

variables  such  that  F.(q.  „  ;f,  )  -  6.  =  o  (1).  Then 
1  i  P, 

VO  -  F1(^,N;r’1)  '  W*”  +  «i  -  . 

The  purpose  of  lemma  4.4  is  to  show  that  assumptions  AIq  and  A2q  imply  a  slightly  more 
general  version  of  assumption  III  used  by  Chibisov  (op.  cit.).  Lemma  4.4  is  then  used  in 
the  proof  of  lemma  4.5  which  is  the  non-null  version  of  lemma  3.2  given  in  Chibisov. 

Lemma  4.4:  Suppose  that  6  and  0  are  any  two  elements  of  0.  For  i  =  let 

_  _  y  _  _ 

S  =  F0(q^(9);9)  and,  for  £  >  0  define  the  sets  B(e)  =  (6 ^-c ,6 ^+e )  and 

S(e)  =  UF'_j{x  :  !fq(x;8)  -  fi.j  <  e}.  If  A10  and  A2g  are  satisfied  then  there 

exists  an  £  >  0  such  that  F^(x;9)  is  differentiable  in  9  at  9  and 

Fq  ( x ;  9  )  =  Ffl(x;9)  +  fQ  ( 9 1  )  (1 ,  x )  ( 9-9  )  +  r(x,0,9) 

where:  i)  v(u)’  =  ( 1 , ( FQ  1  ( u )-8  ] /fi  )f  ( fq 1 ( u )  )  is  bounded  and  continuous 

for  u  in  B(e),  and 

ii)  sup  |r(x,9,9)|  =  o(l9-?Sll)  as  18-91  +  0. 

xeS(c) 

Proof :  Let  0  =  min(61 , 1-0R  .j  )  which  is  positive.  Next,  for  0  <  e  <  0  and  0  <  T  <  ®2 

let  0(e)  =  {x  :  f'V?  -e)  <  ^i+®2x  4  F~’  (S"K_1+e  )} ,  let  M(£)  =  max  (1+x2)/2, 

xeD(e) 

which  is  finite  since  D(e)  is  compact,  let  B($,T)  =  {9  :  I9-Sl  <  t} ,  let 
A  ( £ ,  T  J  =  {y  -  (1,x)0  :  x  is  in  D(e)  and  8  is  in  B(0,T)},  and  let  B(e  ) 
denote  the  closure  of  the  set  B(e)  defined  in  the  statement  of  this  lemma. 

It  follows  from  Alg  and  A2q  that  for  any  e,  with  0  <  E  <  0,  Fq(x;9)  is 
continuously  differentiable  in  9  at  9  for  any  x  in  D(e).  Thus,  by  Taylor's 
theorem  with  Cauchy  remainder,  there  exists  an  a  (0  <  a  <  1),  which  may  depend 
on  x  and  0,  such  that 

F0(x;9)  =  F0(x;9)  +  f Q( ( 1 , x >0 > ( 1 , x ) ( 0-9 )  +  r(x,0,9) 
where  r(x,9,9)  =  ( 1  ,x )  (9-0  )  [f  Q(  ( 1  ,x  )9  )  -  f  ((1,x)9)]  and  0  =  0+  a(9-9). 
Furthermore,  since  0  <  £  <  0  implies  that  B(£)  is  contained  in  the  interval 
(0,1),  and  since,  for  such  e,  v(u)  as  defined  in  (i)  is  continuous,  and  thus 
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(4.21 


(which  equals  O)  and  for  any  0  <  u  <  1, 

F"1<u.-n  =  (K01(u)  -  fi)/f-2 


ar.<i ,  f  or  C  -  0  or  1  , 

.3 


iT  ^  e  O,  J  r,);r ))  ■  -  ♦  r,2gi<OHi,qi<a> 


i  =  1  / 2 , . . . , K- 1 


fgF  F( .(g. (£);■;):  •  *  "‘VW5  1  F  ?2g.  (5))(1,g.<e))  . 

S  i  nc*‘ ,  for  J  =  0  or  1,  f  ^  ( x )  is  continuous  for  any  finite  x  by  assumption  A2^, 
the  next  two  lemmas  are  st ra i ghtf orward  consequences  of  corollary  3  in  Mann  and  Wald  (1943) 
ari  Taylor's  theorem  with  Cauchy  remainder. 

t.emma  I  •  1  :  If  (.  and  ?  are  any  two  elements  of  3,  if  8  -  £  =  0  (1//N),  and  if 

N  P1 

and  A 2?  are  satisfied  for  l  -  0  or  1,  then 

WV’?)  ■  Ft‘9i«e»»?) 

-  (?2A2)fJt(£i  +  ?2g.  <?))<1,g.<£))<9*  -  £) 

+  o  (1//N) 

1 

foi  I.  =  0  or  1  and  i  =  1,2,...,K-1. 


Lemma  4.2:  If  A2n  is  satisfied,  then  F^tx)  is  continuously  differentiable  in  x  and, 
for  0  <  t  <  1  and  0  <  6  <  1, 

F~V)  =  )  +  (t-6  )/f  0(  f"1  (6  ) )  +  r(t,6  ) 

where  sup  |r(t,6)l  =  o(T)  as  t  ♦  0. 

!  t-6  |<x 

The  next  lemma  is  a  basic  fact  involving  the  empirical  process.  The  statement  of  this 
lemma  is  a  slight  modification  of  lemma  3.3  given  in  Chibisov  (op.  cit.)  which  takes  into 
account,  the  fact  that  it  will  be  applied  under  H.j,  not  Hq.  The  proof  of  this  lemma  can 
he  found  in  the  proof  of  Theorem  6  (page  437)  jn  Oihman  and  Skorohod  (1974). 


F.emma  4 . 3 


Suppose  that  F  (xjf.  )  is  continuous  in  x,  that  £  is  an  element  of 


T  -1 


and  let  .  -  F^(q^(^);r.  ).  For  a  >  0  let  B(fl^,a)  = 


v  x  :  I  F 1  f  x  ; 


i  .  1  <  ci}  and  let  'i)  ,c 1)  “  sup  /N|Fk,(x)  -  F,(x?£^)  - 

l  N  l  N  1 


F  (n  (r))  ♦  f  I  .  Then,  for  any  r  >  0 

N*  '  i  i 


xeB(<*>  .  ,(i ) 


1  it  1  i  f*  F  f  <)  (  r  ,  i  )  >  f  }  -  0 . 

i  N  ■ 

t  *0  N 


-  1"- 


Proof  ;  From  (3.6)  and  the  definition  of  6^)^  /£)» 


U<VV  =  ♦  locn'.e^  -  b(2)(e\f>]  -  «°<VV  • 


*  *  — _  m 

Furthermore,  lemma  5.2  shows  that  9  and  9  jointly  satisfy  A  (£,£).  Thus, 

N  N  l 


(5.14)  shows  that 


S°(9*0  )  =  ,$0(5.5)  +  (60-'ft*,O  -  6°<5,S>)  +  [6°(C,0m)  -  «°(5,£))  +  oD  (1//n) 

N  N  N  N 


and,  since  ?2/£2  ~  80  that  n9(?,5)  =  b(p)n°(5,£),  (5.16)  and  (5.17)  yield 


(5.38)  «°<e*,0N)  =  6°(T,o  +  n°(T,c)(i,g(T))t(8N-?)  -  b(p)(e*-oi  +  op  <i//n)  . 


Thus,  since  (5.36)  shows  that 


<0N-f)  -  b(p)(9*-5)  =  (0*-C)  +  E(5  )(n°  )-1  (U(0*,@N)  -  8(2)(c\e)]  +  op  ( 1//N) 


u(0*,qn )  =  61(c,c1)  -  «°(c,l)  +  >  -  b(2) tc1  ,e )i 


-  n°(?,c)(i,g(e))(0*-c>  +  o  d//N)  . 


Therefore,  (5.37)  shows  that 


(5.39)  u(9  .9  )  =  B,,,(5  ,5,1)  +  (M,(C,5)  :  M  (€  .£.?)] 

N  N  (  1 )  7  *8 


F„(g(5))-«1(5,F,1) 


9n  -  ? 


+  o  (1//N) 

1 


and  the  theorem  now  follows  from  corollary  5.1. 


2  •  2  “2 

6.  Finite  Sample  Size  Approximations  to  the  Non-Null  Distributions  of  RJJ ,  R^ ,  and  R^  • 
Each  of  the  statistics  considered  in  this  paper  can  be  expressed  as 

(6.1)  NXjJjB^Xpj 

where  XN  is  a  random  vector  and  is,  in  general,  a  random  positive  definite  matrix. 

The  procedure  employed  here  for  finding  a  finite  sample  size  approximation  to  the  distribu 
tion  of  (6.1)  under  Hi  involves  the  following  steps.  First,  write  as 

(6.2)  Bf,  =  T„T* 


where  Tv  is  a  lower  triangular  matrix.  Then,  find  the  non-null  distribution  of 


Finally,  if  /n  (T^Xn-u)  +  nk-1*°'^'  *  finite  sample  size  approximation  to  the  distribu¬ 

tion  of  (6.1)  is  given  by  the  distribution  of 

(6.4)  Y*  Y 

where  Y  -  Nr_^(/n  u,$)  (see  the  comment  following  theorem  5.1). 

To  accomplish  the  first  step  in  the  procedure  outlined  above,  let  b  *  (b^ , . . . ,b*) 1 
where  bj  >  0  for  j  =  1,2,...,K,  and  define 

(6.5)  Pj(b)  =  1/bj),  j  =  1,2,  ...,K  . 

Some  properties  of  Pj(b)  are  recorded  in  the  following  lemma. 

Lemma  6.1s  i)  Pj(b)  =  bjPj_.,(b)  +  (It^"^)  2  <  j  <  K 

li)  bj  +  bj+1  -  bj[Pj_^(b)/Pj(b)]  -  Pj+1(b)/Pj(b)  2  <  j  <  K-1 

iil)  Pj(b)/Pj..,(b)  =  bj  +  1/b,)-1  2  <  j  <  K 

iv)  Pj.^bJ/Pjtb)  -  ~  -  tb^j^  1/b^'1  2  <  j  <  K  . 

Proof  :  (i)  follows  from  the  fact  that  Pj(b)  =  ( b  j  IT  J=jb^  )  t  ( Vb^  )  +  1/b^)  .  To  show 

(ii)  notice  that  bj  +  bj+1  -  b* [Pj_, (b) /Pj (b) )  »  {bj [Pj (b)-bjPj_1 (b) )  + 
bj  +  1Pj(b)}/Pj(b)  -  {bjfirj'^bj)  +  bj  +  1Pj(b)}/Pj(b)  -  Pj  +  1(b)/P.  (b).  To  see  (iii) 
divide  (i)  by  Pj_-|(b)  and,  finally,  to  obtain  (iv)  divide  (1)  by  P  j  ( b )  and 
rearrange  terms.  ^ 

The  next  lemma  provides  the  principal  result  by  means  of  which  (6.2)  is  accomplished. 

K~  1 

Lemma  6.2:  Define  the  matrix  A  »  <ai,j)i,j=1  by 


a 


ir) 


if 

if 

if 


I  i~1 1  -  0 
I  i~ j I  =  1 
li-jl  >  1 


where  b  »  (b1,...,bK)'  and  bj  >  0  for  j  -  1,...,K.  Then  A  »  TT'  where  T  = 

^i.j’i.j-l  and 
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[P.+1(b)/P  (b)]'*  if  I  i- j |  =  0 


v, 

t.  .  =\  -b. [P.  (b)/P.(b)]/2  if  i-1  =  j  i  =  2,3,...,K-1 

1  1 


if  | i  —  j |  >1  or  j  >  i 


Proof :  Direct  multiplication  of  TT‘  and  part  (ii)  of  lemma  6.1. 

For  0  and  9  in  0  let 


46°(9,9)  = 


<s° ( e ,9 ) 


if  i  =  1 


6°(9,0)  -  (6,9)  if  i  =  2,3, 


1  "  5k-i<0,0) 


if  i  =  K  . 


If,  in  lemma  6.2,  b^  =  1/66 0,0),  then  it  is  shown  in  the  appendix  that  A 
Thus,  it  follows  from  lemma  6.1  that 

(6.7)  D°(0‘,0)"1  =  T{0,0)T(?,0)’ 

where  the  (i,j)th  element  of  T(0,0)  is  given  by 


<1/66°  (?F,0)  +  1/6°(F,8))1/2  |  i-j  |  =  0 


(6.fl)  t.(j<9,0)  =  <  -  ( 1/66?  (0,0)  -  1/6°  (9 ,9  ))1/z  i-1  =  j  i  -  2,3, ...  ,7 


|i-j|  >1  or  j  >  i  , 


"otice  t-hat  when  0=0,  6° (0,9)  =  5°  =  fq(c^).  Therefore,  T(0,9)  =  T  where 


[1/66°  +  ,  +  1/6°]  |  i-j  |  =  0 

[1/66°  -  1/6°]1/2  i-1  =  j  i  =  2,3,...  ,K- 1 


I  i  — j |  >1  or  j  >  i  , 


hd  ’j-'  =  TT’  .  Furthermore,  since  T  ^  =  <  fl  ° )  D  ^(H°), 


(6.10) 


r'1  =  w 


(6.18) 


(6.19) 


a)  y  '  Vl(/iT  v'p(nur'8(2)(5',(i,  v'P(nu)"1t(2)(C1,?)(nu)'1piv) 

if  0  =  8  .  It  should  be  noted  that  when  0  -  8  =  o„  ( 1 ) ,  (5.30)  shows  that 

N  N  N  N 

9)  Ve(6**)  =  T'UfS*,©**)  +  op  (1/Vn)  . 

Thus,  in  this  instance,  the  finite  sample  size  approximations  to  the  non-null  dis- 

tributions  of  *  and  Rj/9N'9N  *  are  the  same-  In  general,  however,  if 

*  ** 

9^  -  0^  /  Op  ( 1 )  then  these  approximations  differ.  In  particular,  lemma  5.2 

shows  that  ^  °p  80  that  the  approximations  to  the  non-null  distribu- 

•  2  ^  * 2  -*  * 

tions  of  and  n'  N*  are  not  the  same# 

Now  consider 


*  **  *  ** 


RK(0N'0N  5  -  NtT(V9N  ),U(0N'0N  ,),(T(9N'0N  ),U(V0N  )] 


and  notice  that 


T(V9r,,u(<VeN*)  "  T(c,c)*B(1)(c1,e,l) 


(6.21)  T  ( 0  *  ,  0*  *  )'  [U(9*,9**)  -  S(  ,,  (C1  ,C,C  )] 


+  tT(0*,0**),6(1)(C1,C,!)  -  T(C,?)*0(1)(S1,C,5)1  . 

***  *  ** 

Thus,  the  limit  distribution  of  t^®n'9N  *  under  H1  does  not  follow  directly 

from  theorem  5.1. 

•  **  0 

Lemma  6,3:  Let  9^  and  9^  denote  two  estimates  of  6  and  suppose  that  there  exist 

—  ~  *  *  *  _ 

5  and  5  in  H  such  that  0N  and  0^  jointly  satisfy  A^(5,5)>  Further¬ 

more,  let  J ( 8 , 9 )  denote  the  Jacobian  of  the  transformation  from  6° (0,8)  to 
<*(9,9)  =  T(  9 ,0  ) 1  1  j  (5  ,5,5)-  If  Alj  and  A2^  are  satisfied  for  l  “  0,1 
then,  under  H^, 

t(9*,9**),B(„(51,T,C)  -  T(I,C)'0{1)(C1,C,l)  +  J<5,5)  [6°(e*,0**)  -  «°(T,C)3  +  op  (1//N)  . 

0***  n  —  —  . — 

Proof :  Under  the  hypotheses  of  this  lemma,  6  B  (5,5)  «  0p  (1//N)  (see  the 

proof  of  theorem  5.1).  Furthermore,  from  the  exact  expression  for  J(8,9)  given 


in  the  appendix,  it  is  seen  that  its  elements  are  continuous  in  6  (8,9).  Thus, 


the  lemma  follows  by  the  standard  arguments  due  to  Mann  and  Wald  (op,  cit.). 


Theorem  6.1;  Let  8^  and  8^  denote  two  estimates  of  6'  and  suppose  that  there  exist 

_  *  «•  —  - 
£  and  £  in  3  such  that  8N  and  8N  jointly  satisfy  A^  (£,£)•  Define 

B  (£\f.?>  as  in  Theorem  5.1  and  let 

Mg(£1,r,S)  t(T(£,£)>  -  J(£,£))n°(£,£)  -  T(£,i)'iT1<£#£1)ni,g<£)J  . 

«1B(U)  -  -(T(£,£)'  -  J(?,£)]n°(£,£)  [1,g(£>]  , 


$(5)<c\s.?>  -  [T(£,£>*  :  Mg<£\£,£)  :  M10(£,£)]S1(£1,£/£) 

x  [T(T<C  ) '  •  Mg  (£ 1  ,T»5 )  ;  M1(j (£",£)] '  . 

If  Alj  and  A2g  are  satisfied  for  l  »  0,1  then,  under  H^, 

/N[T(0*,0**)'U(8*,0**)  -  T(T,?),B(1)(c\T.?)]  +  Nx_i(°4(5)(£1J,£))  . 

Proof i  It  follows  from  lemma  6.3,  (6.21),  (5.14),  (5.16),  (5.17)  and  (5.18)  that 

T(8*,9*Vu<e*,e**)  -  T(T,l)*e{1)(C1,T,l)  -  t(?,£>* [u(8*,e**>  -  e(1)(£\F,?>] 

+  J(£,£)[S°(8*,e**)  -  «°(£,£)]  +  o^  (1//H)  - 

t(T,|)*[i;  M^c’.e.c) :  M2(£,£)]((FN(g(7)>  -  s\T  c1))*  \  ( e*-T) •  •  <e**-c > * r 

♦  J(£,£)[-n°(£,£)[1,g<e)]  i  n°(£,£)(1,g(T)]]  [(6*-0’  (9**-? ) ' )  *  +  op  (1//N)  . 

Since  Mgtc1,!,!)  -  T(?.£)’M1<£1,?,£>  -  J (£,£  )n° ( £  ,£  )  [1,g<£)]  and 

T(£,£)’M2<£,£)  +  J(T.£)n0(£,£Hl,g<£>], 

T(8*,8**),0(8*,e**)  -  T(f,£)'B(1)<£ \?,£)  - 

[T(£,£>*  :  Mg(£\l,£)  :  M10(U)1  I(FN(g(?))  -  6  1  (T,C 1 )  >  ’  ;  (9*~T>  *  ;  (6**-£)']  '  +  op  (1//n) 

and  the  theorem  follows  from  lemma  5.1.  ^ 

•  0  — 
Theorem  6.2 ;  Let  8n  denote  an  estimate  of  8  and  suppose  that  there  exists  a  £  in 

*  —  » 

=  such  that  0^  satisfies  condition  Define  £  as  in  lemma  5.2  and 

define  B^1j(£1,£,f)  as  in  theorem  5.1.  Let 
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=T(£,£)'  -  [T(£,£)'  -  J(C,C)]n  <£,£)  [l,g(5)]E(Otn  ) 

=  -[M11(T.?»n1(T,c1>  +  (Ttc.e)’  -  j<5,i)>n°(c,f:)]<i,g<f)]  , 


If  A1  ^  and  A2^  are  satisfied  for  £  =  0,1  and  if  0^  is  defined  by  (3.15) 
then,  under  Hf# 

^[T(VV'I,IVV  ‘  T<^),8(1)<s1,£.!)]  *  6)u',T'F,))  • 

Proof :  It  follows  from  (5.36),  (5.37),  and  (5.38)  that  ~  6^(£,£)  ”  Op  (1//N)  + 

ir°(T,I)[i,g(I)][E(T)(n0)"1  |  i-E(t)(n°)‘1ii1(T,e1) M,g(?>) J [<FN<g<f)>  -  6 1  1 ) > • 

•(e*-?)1]’  • 

Thus,  the  theorem  follows,  after  some  rearranging  of  terms,  from  (5.39)  and 

corollary  5.1  as  in  the  proof  of  theorem  6.2.  ■ 

2  *  •* 

The  finite  sample  size  approximations  to  the  non-null  distributions  of  ) 

ard  R^(0n,S)  follow  from  theorems  6.1  and  6.2.  They  are  given  by  the  distribution  of 


(6.4)  where 


Y  ~  Vl(/N  t^.O’B(1)(£1,?,C),  t(3)(£1,r,?)) 


for  general  choices  of  0N  ,  and 


y  -  Vl(/N  T(£,£)I8(1)(£1,r,£), 


wh.en  =  9 

N  N 


The  following  lemma,  which  is  easily  demonstrated  and  covers  all  of  the  cases  in  this 
section,  provides  a  convenient  expression  for  the  distribution  of  (6.4). 

Lemma  6.4:  Suppose  that  Y  ~  N  (w,$),  the  the  rank  of  t  is  r  4  p,  and  that  u  is  an 

element  of  the  space  spanned  by  the  columns  of  $.  Let  >  X^  >...>  X  >  0 

denote  the  non-zero  characteristic  values  of  $,  let  A  =  diag(X ^ , . . . ,X^ ) , 

let  P  denote  a  p  «  r  matrix  which  satisfies  $P  ~  PA  and  P'P  =  I, 

J/, 

and  let  w  =  A  ^P'u.  Then  the  distribution  of  Y'Y  is  the  same  as  the 
distribution  of 


ri=iVW‘ 


where  Z.| ,  Z2  , .  . .  ,Zr  are  independent  standard  normal  random  variables. 

By  means  of  Laguerre  series  expansions  (see  Gurland  (1953,  1955,  1956),  Kotz  et  al  (1967), 
and  Gideon  and  Gurland  (1976a,  1976b)),  percentiles  for  the  distribution  of  (t>.23)  are 
easily  computed,  and  finite  sample  size  approximations  to  the  power  of  the  tests  based  on 
the  statistics  considered  here  are  easily  obtained. 

7.  An  Extension. 

While  the  results  of  the  previous  sections  can  be  generalized  in  many  ways  to 
accomodate  other  types  of  parameters,  it  is  interesting  to  note  that  the  following 
extension,  which  includes  several  useful  examples,  involves  only  minor  changes  of  notation 
in  the  previous  arguments.  In  place  of  A1^  and  A2^  (4  =  0,1),  assume: 

I:  There  exists  a  parametrization  of  F^(x)0)  so  that 

i)  0  -  (0  -  <VV  :  <  01  <  02  >  0} 

ii)  Fq(x;0)  =  Qq(0  1 ,02<tQ(x) )  where  Qg  maps  R  onto  the  interval  [0,1]  and 
is  a  horaeomorphism  mapping  a  set  Ag  (which  is  free  of  0)  onto  R. 

II:  Qg(y)  is  continuously  differentiable  and  gD(y)  *  ^-Og(y)  is  non_zero  f°r 

finite  y. 

Ills  There  exists  a  parametrization  of  F^(xsC)  so  that 

i)  =  -  u  -  <e , ,c2 > •  :  —  <  e,  <  ?2  >  0} 

ii)  F  ^  ( x  j  C )  “  +  Cj't'^x))  where  Q1  maps  R  onto  the  interval  [0,1]  and 

4^  maps  a  set  A^ ,  which  is  free  of  £  and  contains  A0,  onto  a  set  B1 . 

IV:  Q 1 ( y )  is  continuously  differentiable  and  qj(y)  “  njy  P 1 ( y )  is  non-zero  for  all 

finite  y. 

V:  For  l  =  0,1,  $^(x)  is  continuously  differentiable  for  all  x  in  Ag  and 

T^(x)  =  <^(x)  is  non-zero  for  all  x  in  Ag . 

It  should  be  noted  that  the  assumptions  that  qj  and  x1  be  non-zero  are,  in  general,  not 
necessary.  They  are  included  here  for  convenience  only  to  ensure  that  the  ranks  of  the 
matrices  i  =  1,...,6,  remain  unchanged. 

The  following  examples  satisfy  assumptions  I,  II,  and  V. 


Example  1 : 

(7.1) 

I 

Example  2; 

(7.2) 

Furthermore,  if  one  of  these  examples  'ere  chosen  to  represent  Hq,  then  since  the  other 
example,  as  well  as  any  F,  satisfying  Al,  and  A2,,  satisfy  III,  IV,  and  V,  these  would 
be  candidates  for  an  alternative  hypotheses  in  the  type  of  comparison  considered  here. 

Assumptions  I,  II,  III,  IV,  and  V  require  the  following  changes  in  section  3  through 
6.  First,  the  definition  of  ( 0 )  (see  the  comment  following  (3.3))  becomes 


(7.3) 

g^g)  =  *Q1(vi(0))  i  *  1,...,K-1  , 

where  v . (9 ) 

1 

=  (c. -9,)/92  and  -  Q01<5i)"  Next'  let 

Y± ( e  )  -  $, (♦g1 (vi (5 )  ) )  i  =  1,2,...,K-1 

(7.4) 

and  let 

(7.5) 

ni<5)  =  T1<*01<Vi(5,))/T0('>01(vi(5)))  1  = 

Then,  (5.3), 

(5.4),  (5.7),  (5.8),  (5.11),  and  (5.12)  become 

«°("9,9)  =  90(91+92vi(?)) 

(T,C )  -  0,(5,  +  52y.(5)> 

"°(9,0>  =  q0(91  +  9^(9)) 

o  o  —  — 

(7.6  ) 

",  -  ",(9,9)  -  q0<ci) 

nj(5,5)  =  q  1  (5  ,  +  52Y1(E))hi(5) 

|a  F  (g  (9);9)  =  K° (?  ,0  )  ( 1 , v  ( ¥ ) )  • 

a  n  U  1  1  1 

Fo(gi(n);0)  =  (-02/e2)r°(9,9)(1,vi(e))'  =  -*  °  (9  ,0)  (  1 ,  v  .  ( 9  )  )  1 


The  log-normal  distribution  with  density 

p(x;u,o)  =  [/2x  o)  1  |  exp{-  |  ( *  n-q*  — ) Z}  x  >  0 
where  -*•  <  u  <  +°°  and  a  >  Q.  Here  ( y  >  =  ${y),  where  ♦(•  )  denotes  the 
distribution  function  for  a  standard  normal  random  variable,  <i0(x)  ”  In  *» 

6  ,  =  -u/0,  8^  =  1/0,  and  A^  =  (x:x>0). 

The  Weibull  distribution  with  density 

Yl'1  Y1 

P(x|yi(Y2)  -  (Y^YjHx/Yj)  exp{-(x/Y2>  }  x  >  0 

where  Y1  and  y2  >  0-  Here  QQ(y)  =  1  ~  exp{-ey} ,  ♦gf’O  *  In  x, 

9 1  =  -Y ,ln  Y2,  9 2  =  Y  , ,  and  AQ  =  {x  :  x  >  0). 
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and 


(7.7)  It  F,(g.  (5)sU  -  <-C,/C,)*  V,.F,)<i,v1(0)‘  =  -!'(5,5)(1,v,(E)]'  . 

d  S  1  1  ^  jC  1  1  1  1 

Finally,  except  where  the  function  g^  plays  the  role  of  an  argument  of  FN(x),  replace 
g^  by  v^ .  Then  the  arguments  of  the  previous  sections  are  unchanged  and  all  results 


given  there  follow  from  assumptions  I  through  V. 


The  two  points  to  be  covered  are:  i)  exhibit  the  matrix  J(8,0)  for  0  and  0  in 
0,  and  ii)  show  that  the  quadratic  form  (3.8)  is  equal  to  the  standard  definition  of  the 


Pearson  chi-square  test  statistic  when  unknown  parameters  are  present  and  the  interval 
boundary  points  depend  on  the  unknown  value  of  0. 

Let  *  *  ($  1 , . . .  )  *  and,  for  0  in  0,  let  io(0,0)  =  T(0,9)'<t>  where  T  ( 0  , 9  ) 

is  given  by  (6.8).  Then,  recalling  the  notation  in  (6.6), 


n/A6°  +  1(¥,0)  +  1/6°  {0 ,9  )  ]V2  <»>i 


U5±  (0,0) 


(1/A6°+1(?,0)  -  1/«°+1(0,0)]1/2$.  +  1  i  =  1 » 2  , , . .  ,  K-  2 


( 1/A6°  (0 ,0  )  +  1/«°_1(0,0)]1/2I*’k-1  i  =  K-1  . 


_  (0,0)  K-1 

Thus,  J(0,9)  =  ( — tT - )  and,  if 

3«“(0,0)  i,j“1 


ri<0,9)  =  1/A6°+1(?,0)  +  1/6° (0 ,0 ) 


r^ (8,8  J  -  1/A6°+1 (9,8)  -  1/d°(9,0) 


<  0 ,9 )  =  1/A«°+1(0,0)  -  1/6°+1(8,9) 


Si(0,9)  =  (1/A«°  +  1(9,9))2  -  ( 1/6° (0 ,0  )  )2 


(0,8)  =  ( 1/A6  °+1 (0 , 0 ) ) 2  -  (1/i°+1(8,0))2 


S^F,©)  -  (1/A6°+1(8,8))2 


lenient  of  J(8,9)  is  given  by 


-l/o  -  -  -  -  Vo'  - 

i  48  !  2  S  (8 ,0  )*  -r.  (9,8  )  /2  S,  (8 


S.  (8 ,0 )$ -r.  (9,8  )  Si (8,8  )*i  +  1)  | i- j [  =0  i=1 ,2, . . . ,K-2 


i  =  j  =  K-1 


•1  'r  1^,8)'  1/2  S.  9  >$.-r.  (8,8  )“  1/2  S.  (0,9)4.  ]  i+1  =  j  i  =  1/...,K-2 

l  l  ii  l  1  +  1 


in  I i-1 |  >  1  or  i  <  j 

1  the  definition  of  Y,(9,8)  given  in  (5.24).  Then,  the  standard  definition  of 
iv  ii  son  -hi-square  test  statistic  is  given  in  terms  of  the  quadratic  form 

X^(8,8)  =  n  E^=1Y. (8,8)2/A6°(8,8) 

A  .  5  ' 

—  •  —  -1  — 

=  N  Y(0,0)*D(0,0)  Y(  0 ,9 ) 

_  —  —  .  _  0—  0  — 

whore  Y'0,0)  =  (  Y1  (  9 , 9  )  , . .  .  ,  Y^Q  ,  0  )  )  1 ,  D(0,0)  -  diag(AA  1  (0 ,0  > , . . .  ,A6R(0  ,0  ) ,  and  0  and 

“  are  elements  of  0- 

To  show  that  (A.5)  =  (3.8)  let  AK  denote  the  K  x  K  matrix  with  ones  on  and  below 

_  Q  _  ]£_  J 

the  main  diagonal  and  zeros  elsewhere.  Also,  let  G(9,0)  —  ^  t  j  )  (®  f  j-j  an(^  no^e 

that  n°(9,9)  =  G(8",8)  -  6°  (IT  ,8  )6°  ("ft",  8 )  1  (see  (3.7)).  Then,  it  follows  from  (5.24)  that 
AyYtO.fl)  =  (0(0,8)  • ,0] *  and,  therefore, 

(A.fi  X^(?,8)  =  (UlfMr.OHA^te^JA^)  [udMr.O]’  . 


(AKD<9,9 )A^) 


G(9 ,0  )  6° ( 8 ,8  ) 


5°  (8 ,8)  1 


(G-fiV')"1 


-  ( G-5  °5  0 1 ) ” 1 6  0 


(A  D( 9 , 8 ) A ' ) 
K  K 


,0,0.  -1  1+6°' (G-6°50' )-1<5° 


where,  for  convenience,  the  dependence  on  8  and  8  has  been  surpressed,  and  X^(0,0)  ** 
R^(?,9 )  hy  (A. ft)  and  (A. 7)  . 

A  convenient  expression  for  D^(0,0)  1  is  obtained  by  noting  that  {A^D(9,9)A^)  = 

( A^)“1f‘(0^#0  j"1  ( Ax  )~1 .  Since  (Aj^)”1  is  the  K  x  K  matrix  with  ones  on  the  main  diagonal 
and  negative  ones  on  the  diagonal  just  below  the  main  diagonal,  and  since  D(0,0)  = 

diagl  1/A<5°('0,e),...,1/A«°(0<0))/  (  A^f  0  ,9  )A^  )-1  =  <qi(.>*j=1  where 


Thus , 
(A. 8) 


D°(6,9)'1 


"i.S 


K-1 

i,j=1  ’ 
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